ON THE SEMI-CLASSICAL LIMIT FOR THE NONLINEAR 
SCHRODINGER EQUATION 

REMI CARLES 

Abstract. We review some results concerning the semi-classical limit for the 
nonlinear Schrodinger equation, with or without an external potential. We 
consider initial data which are either of the WKB type, or very concentrated 
as the semi-classical parameter goes to zero. We sketch the techniques used 
according to various frameworks, and point out some open problems. 
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1. Introduction 



Consider the nonlinear Schrodinger equation (NLS): 



(1.1) 




where the potential V = V{x) and the nonlinearity / are real-valued. In some spec- 
ified cases, we allow the potential to be time-dependent. To simplify the discussion, 
we assume that n ^ is an integer. More precise assumptions will be made accord- 
ing to the different cases we study. We assume e e]0, 1], and we aim at describing 
the asymptotic behavior of u £ as e — > 0, for the following two families of initial data: 

Monokinetic WKB initial data: 



(1.2) u £ (Q,x) = a e {x)e iMx)/£ , with a £ (x) ~ a (x) + ea x {x) + e 2 a 2 {x) + . . . , 



for some point (xq,£o) m the phase space M 2 ™, independent of e. 

There are at least two motivations for such a study, referred to as semi-classical 
analysis or geometrical optics. We outline them here, and refer to the survey [47] 
for a broader discussion on this subject. The first one comes from the applied 
mathematics, and may find its origins in physics. In the case of suppose that 

e represents the (rescaled) Planck constant. It may be small compared to the other 
parameters at stake. In this case, it is sensible to consider that the asymptotic 
behavior of u e as e — > provides a reliable approximation of the exact solution. 
Hopefully, the asymptotic model is easier to describe than the initial one (jl.lj) 
(|1.2p . If V is a confining potential (e.g. harmonic potential), then (jl.ip may be a 
model to describe Bose-Einstein condensation; see for instance [22l [46] . The value 
of k then depends on the asymptotic regime considered. Another motivation stems 
from the propagation of singularities for equations where the small parameter e 
is not necessarily present initially. Most of the studies in this direction concern 
hyperbolic equations. However, this field is applicable to Schrodinger equations as 
well (see e.g. [TJ [4TJ [51] ) . The following illustration is a straightforward consequence 
of the analysis presented in §3.21 

Theorem 1.1 ([13], Cor. 1.7). Let n > 3. Consider the cubic, defocusing NLS: 

(1.4) id t u + ^Au = \u\ 2 u, x <E 1R" ; u\ t=Q = u . 

Denote s c = § — 1. Let < s < s c . We can find a family (uq)o< £ ^i in <S(R™) with 

II«oIIh s (K") ^0 as e -> 0, 
and < t e — > such that the solution u £ to (| 1 .4[) associated to Uq satisfies: 



in the sense of asymptotic expansion. 



Concentrated initial data: 



(1.3) 




\\ u£ (t £ )\\H k {R™) — * +°° as £ — > , Vfc G 



s 



J 2 
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This result was first established in [21] in the case k — s. The fact that one can 
consider a broader range for k, in the spirit of [42], relies on a fine analysis of the 
limit for (|l.ip - (|1.2|) , provided essentially in [52"] . 

1.1. Monokinetic WKB initial data. In the case of initial data of the form 
(|1.2p , an approximation of the form 

(1.5) u £ (t,x) ~ (ao(t,a:) + ea.i(t t x) + s 2 a 2 (t,x) + . . .) j**.*'*)/* 

e— >0 

is expected. Note that only one phase and one harmonic are sought: this is an 
important feature of Schrodinger equations with gauge invariant nonlinearity. In 
the case of wave equations for instance, the story is completely different (see e.g. [J7J 
and references therein). Note also that such an approximation must be expected for 
bounded time only. Even in the linear case / = 0, a caustic appears in finite time 
in general. Near a caustic, all the terms ao, ai, ...become singular. Past the 
caustic, several phases are necessary in general to describe the asymptotic behavior 
of the solution (see e.g. 25J for a general theory in the linear case). However, we 
will see that the analogous phenomenon in the nonlinear setting (say, f(y) = y) 
with k = (highly nonlinear regime) might be very different. 

Plug a formal expansion of the form (|1.5p into (|1.1[) . Ordering the terms in 
powers of e, and canceling the cascade of equations thus obtained is aimed at 
yielding a , a 1; . . . 

Assume for a while that k ^ 1. To cancel the term of order 0(e ), we find 
ao (d t $ + i | V$| 2 + Vj = ; $, t=0 = </>o • 

Since we seek a non-trivial profile ao, we impose a stronger condition: $ must solve 
the eikonal equation 

9t $ + I| V $| 2 + y = ; $| i=o = 0o- 

Canceling the term of order 0(e 1 ), we get: 

1 { if k > 1, 
5 t a + V$ • Va + -a A$ = { . . „ 

2 [ - 2 /(|a | 2 )a j£k = 1. 

We see that the value k = 1 is critical as far as nonlinear effects are concerned: if 
k > 1, no nonlinear effect is expected at leading order, since formally, u £ ~ aoe'*/ £ , 
where $ and ao do not depend on the nonlinearity /. If K = 1, then ao solves a 
nonlinear equation involving /. 

We will see in Section [2] that when k = 1 , ao solves a transport equation that 
turns out to be a ordinary differential equation along the rays of geometrical op- 
tics, as is usual in the hyperbolic case (see e.g. [47]) • More typical of Schrodinger 
equation is the fact that this ordinary differential equation can be solved explicitly. 

Assume now k = 0, and proceed the same way. Plugging (|1.5p into (|l.ip . we get: 

ll — ' " /(|a | 2 )=0, 

V$ • Va + -a A$ = -2if (|a | 2 ) Re (a aT) a . 

We see that there is a strong coupling between the phase and the main amplitude: 
a is present in the equation for $. In addition, the above system is not closed: $ 



(1.6) 
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is determined in function of ao, and ao is determined in function of ai. Even if we 
pursued the cascade of equations, this phenomenon would remain: no matter how 
many terms are computed, the system is never closed (see 30 ). This is a typical 
feature of supercritical cases in nonlinear geometrical optics (see [19l[20]). We shall 
call the study of this case highly nonlinear WKB analysis. We will see in Sj3j some 
ways to overcome the difficulties pointed out above, especially in the case /' > 
(defocusing, cubic at the origin, nonlinearity) . 

Remark 1.2. We consider only monokinetic initial data. Studying the nonlinear 
effects relevant at leading order (re — or 1) when the datum is of the form 

u e (0, x) = a (x)e^/ £ + b (x)e^ x ^ £ (cj> ± po), 

for instance, seems to be an open problem. 



1.2. Concentrated initial data. For data of the form (|1.3p . a formal analysis 
shows that the case re = is critical: if re > (not necessarily an integer), no non- 
linear effect is expected at leading order. We shall therefore restrict our attention 
to the case re = 0. We also consider the case of a pure power nonlinearity, 

/(KI 2 )=aki 2ct , 

for some a > and A S R. In this case, setting u e — e~ n / 2 u e ', (| 1 . lj) — (| 1 .3|) is 
equivalent to: 

(1.7) ie9 t u E + yAu E = W + A£ nCT |u e | 2(T u e ; u e (0,x) = -^R e lx<o/e ■ 

In the case A > (defocusing nonlinearity), dispersive effects are expected to alter 
the concentrated form of the initial data. This is proved in [TTJES] when the external 
potential is a polynomial of degree at most two. It seems that proving a similar 
result in the more general (and fairly natural) framework of smooth, sub-quadratic, 
potentials, is still an open problem. Note also that the dispersive effect can be just 
the first step of the dynamics. It can be followed by a linear dynamics induced by 
the potential. In this regime, the potential may cause a refocusing phenomenon. 
This is the case for instance when V is an isotropic harmonic potential [11J . We 
discuss more precisely these results in £ )4.H 

When A < (focusing nonlinearity), several papers have considered the case 
when the profile R is the ground state associated to NLS without potential, that is 
when R — Q, where Q is the unique positive, radially symmetric, solution of: 

-^AQ + Q + \\Q\ 2a Q = 0. 

When a < 2/n (sub-critical case at the L 2 level) , orbital stability of the solitary 
wave suggests that the solution if evolves under the form 

u ^ t ,x) = Q {^Y^j e fa -««/V^W. 

We will see that this is the case, with (x(t),£(t)) given by the Hamiltonian flow 
associated to — + V: the additional purely time dependent phase shift ip e is 
known explicitly in the case without potential, but not in general. The first mathe- 
matical result on this problem is due to J. Bronski and R. Jerrard [B]. Refinements 
were then given by S. Keraani [35J ISZl I3H] ■ We outline the approach of [35] in ^4.21 
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Note also that the semi-classical limit e — ► for (|1 .7|) is analogous to the long time 
behavior for the solutions to (|1.7j) with e = 1; see e.g. [27] 154"]. 

2. WKB ANALYSIS FOR A WEAK NONLINEARITY 

When k ^ 1, the first step in the WKB analysis presented in HI. II consists in 
solving the eikonal equation. This step relies on the Hamilton- Jacobi theory. It 
is well-known, at least when the potential V and the initial phase 4>q are smooth, 
that the local inversion theorem yields a local in time, smooth solution in the 
neighborhood of (t = 0,x), for all x € R™ (see e.g. [23j). In order to have a local 
existence time which is uniform with respect to x € R™, the following assumption 
is essentially necessary (see e.g. [14]): 

Assumption 2.1. The potential V may depend on time: V = V{t, x). We assume 
that the potential and the initial phase are smooth and sub- quadratic: 

• V e C°°(R t x R"), and d*V € C(R t ; L°°(R")) as soon as \a\ ^ 2. 

• 4> e C°°(R n ), and <9£</> € L°°(R n ) as soon as \a\ ^ 2. 

Remark 2.2. Of course, if we worked on a compact set instead of R™, the above 
assumptions would not be necessary. 

A global inversion result (see [49] or [23]) and Gronwall lemma yield: 

Lemma 2.3. Under Assumption \2.1[ there exist T > and a unique solution 
</>oik € C°°([0, T] x K n ) to: 

(2.1) cWeik + ^|V0 oik | 2 + V = ; <£ oik | t= o = 0o ■ 

This solution is subquadratic: <9"0oik £ i°°([0, T] x R") as soon as \a\ 2. 

Remark 2.4. In [14], examples are given, that show that if either the potential V 
or the initial phase 4>q has a super-quadratic growth at infinity, the above result 
fails. Sub-quadratic potentials play a special role in the mathematical analysis of 
Schrodinger equations: the results of [2"51l2"9"| imply local in time Strichartz estimates 
for the semi-group associated to —A + V. On the other hand, in space dimension 
n = 1, — d% — x 4 is not essentially self-adjoint on Co°(R) (see [26] Chap. 13, Sect. 6, 
Cor. 22]). If V tends to +oo at infinity, with super-quadratic growth, the available 
results are very different from those of the sub-quadratic case, see e.g. [551 155] . 

To prove this lemma, we introduce the Hamiltonian flow: 

f d t x(t,y) =£(t,y) ; x(0,y)=y, 

\dt€(t,y) = -VV{t t x(t,v)) ; f(0,iO = Wo(tf). 

The time T is such that the map y i— » x(t, y) is a diffeomorphism of R" for t £ [0, T]. 
Therefore, the Jacobi determinant 

J t (y) = deW y x(t,y), 

is bounded from above, and from below away from zero, for t € [0,T]. The justifi- 
cation of the leading order asymptotics sketched in ijl.ll is: 

Proposition 2.5. Let k ^ 1 and f € C°°(R+;R). Assume that there exists a 
smooth function ao independent of e such that 

a £ -> a m H s {R n ), Vs ^ 0. 
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Then under Assumption \2.ll for all e g]0, 1], Ijl.ljl — <j 1 - 2 j> has a unique solution 
u e £ C°°([0,T}xR n )nC{[0,T};H s ) for alls > n/2, where T is given by Lemma\EM 
Moreover, there exist a, G € C°°([0, T] x R"), independent of e G]0, 1], where a S 
C([0,T];Z 2 n L°°), and G is real-valued with G € C([0, T];L°°), such that: 



' L°°([0,T]:L 2 nL°°) 

The profile a solves the initial value problem: 



as £ 



(2.3) 

and G depends nonlinearly on a 
a(t, x) 



d t a + V</»cik • Va + -aA^ e ik = ; a\ t =o = oq, 



1 



y/ J Mt,x)) 



ao (y (*,»)) 



G{t,x) = - { /(J.d/^xjj-^aod/ft.a;))! 2 ) 



In particular, if k > 1, then 

1 - nJteik/e 



as e — > 0, 



L°°([0,T];L 2 nL°°) 

and no nonlinear effect is present in the leading order behavior of u £ . If K = 1, 
nonlinear effects are present at leading order, measured by G. 

We see that the critical nonlinear effect (case k = 1) is a self-modulation of the 
amplitude. In the context of laser physics, this phenomenon is known as phase 
self-modulation (see e.g. [57l l4"l [24] ). 

Sketch of the proof. The proof given in [14j consists in changing the unknown func- 
tion, by setting 

a e = u £ e-^ /e , 



where </> c i k is given by Lemma l2~3l Then (|l.ip ~ ()1.2|) is equivalent to: 

7(" -)« ; 



d t a e + V^cik ■ Va £ + ^a e A0 oik = i^Aa e - / .:' 



au_ 



{=0 



Energy estimates show that the above equation has a unique, smooth solution 
a e £ C([0,T];H s ) for all s > n/2, uniformly bounded for e e]0, 1]. This step uses 
the facts that ^ e jk is sub-quadratic and iA is skew-symmetric. We can then neglect 
the terms eAa £ and a§ — clq, so that \\a £ — a e ||z,<*>([o,T];H s ) = where: 



(2.4) d t a e + V0 oik ■ Va e + ia e A0 cik = -fe" 
Recall that Jt{y) is the Jacobi determinant. Denote 



7(|s £ l 2 )s £ 



l |t=0 ^ fl o- 



^ £ (i,j/) :=a^t,x(t,y))^/Myj. 

We see that so long as y i— » :r(t, y) defines a global diffeomorphism (which is guar- 
anteed for i e [0, T] by construction), (|2.4[) is equivalent to: 



L / (j t {y)- 1 \A e \ 2 ) A e ; A e (0,y) = a o (j/). 
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This ordinary differential equation along the rays of geometrical optics can be solved 
explicitly, after we have remarked the identity <9t|A e | 2 = 0: 

A £ (t,y) = a (y)cxp (-ie^ 1 J f (j^ 1 \a Q {y)f) ds^j . 
Back to the initial solution u £ , this yields the proposition. □ 
Remark 2.6. A similar result is proved in [15j for the equation 

2 

ied t u £ + £ —Au £ = V{x)u £ + Vr (J) u £ + \e\u £ \ 2a u £ , 

where Vr is lattice-periodic. The presence of this rapidly oscillatory potential 
changes dramatically the geometry of the propagation. Using the corresponding 
Bloch theory, a similar phase self-modulation phenomenon is proved, under the 
assumption that the initial data are well-prepared. Removing this assumption, or 
considering highly nonlinear regimes (as in @ & re interesting open questions, and 
have physical motivations in the context of Bose-Einstein condensation. 



3. Highly nonlinear WKB analysis: k = 

We saw in Hl.ll that constructing a formal asymptotic expansion for (| 1 . (|1 .2|) 
is a delicate issue when n = 0. We also point out that another problem arises, even 
if one has managed to construct an approximate solution v £ that solves 

(3.1) ied t v e + jAv e = Vv e + f(\v £ \ 2 )v e +e N r s N ; vf t=Q = wf t=0 , 

where N is large, and r e N is bounded in L 2 for instance. Setting w £ — u £ — v e , and 
supposing that u £ and v £ remain bounded in L°°(M. n ) on a time interval [0,t], the 
usual L 2 estimate for Schrodinger equations yields: 

s\\w £ (t)\\ L2 < C [ \\w £ (T)\\ L ,dT + 2e N f \\r £ N {T)\\ L ,dT. 
Jo Jo 

We infer, using Gronwall lemma: 

\\w £ (t)\\ L 2 sC Ce N - 1 e ctl£ . 

The exponential factor shows that this method may yield interesting results only 
up to time of the order ce\ loge| e for some c, 9 > 0. Note that in some functional 
analysis contexts, this may be satisfactory (see e.g. [21], or the appendices in [8lfl3]). 
However, it seems reasonable to wish to have a description of the solution of (|1.1|) 
(|1.2p at least on a time interval independent of e. We list below several approaches 
that yield such information. 

Remark 3.1. In a slightly different context, a fairly explicit example in [12] shows 
that one may find a function satisfying (|3.ip for N arbitrarily large, such that 

liminf \\u £ (t £ ) - v £ (t £ )\\ L 2 > 0, for t £ = e and < < 1. 

E— >0 



Therefore, the stability issue in this highly nonlinear regime is really delicate. 
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3.1. Modulated energy functional. A general technique was introduced by 
Y. Brenier in [5] . It yields the convergence of some physically important quantities 
(such as the Wigner measure, see e.g. [3TJIII]), but not of the wave function u e it- 
self. In the case of the nonlinear Schrodinger equation, it has been used by P. Zhang 
[58] (see also [59] for the case of the Schrodinger-Poisson equation). More recently 
F. Lin and P. Zhang have adapted this approach in the case of the Gross-Pitaevskii 
equation, in the exterior of an obstacle [43] , We shall present the technique of 
Brenier in the case of l|l.l[) - (|1.2[) . using the simplified approach of [35]. In all this 
paragraph, we will assume V = 0: no external potential is present. 

The first step consists in guessing a suitable approximate solution. Even though 
the system (|1.6[) is not closed, the analysis of fJ2] shows that so long as $ is smooth 
and V$ is a global diffeomorphism, the second equation of (jl.6[) is of the form: 

a = iSa , 

where ao stands for the differentiation along the rays associated to V$, and H is 
real-valued. In particular, the modulus of ao is constant along these rays. Setting 
(p,v) = ( | a.p | 2 , V$) as a new unknown function, (|1.6|) yields: 



(3.2) 



d t v + v ■ Vv + f'(p)Vp = ; W| t=o = V0 o , 
dtp + v ■ Vp + pdivv = ; ^^=0 = |o- | 2 - 



If /' > 0, we get a compressible Euler equation, which is hyperbolic symmetric in 
the sense of Friedrichs. We shall assume now that /' = 1, that is, we consider a 
cubic, defocusing nonlinearity in Note that older formal approaches suggest 

the introduction of (|3.2[) as a limiting equation. In [301 Chap. Ill], we find: 



(3.3) 



3 f <F + i|V<F| 2 + |a £ | 2 = £ 2 ^ ; $f t=Q = 
d t a. £ + V$ e ■ Va £ + i a e A$ £ = ; a[ t=0 = a 



Of course, this choice is not adapted when the amplitude a e vanishes, so it must be 
left out for a rigorous mathematical analysis, when G _L 2 (R"). Passing formally 
to the limit e — > 0, the right hand side of the equation for <i> e vanishes, and using 
the hydrodynamical variables as above, we retrieve 



The modulated energy functional associated to (|1.1|) - (|1.2|) when V = and 
f(y) = V is: 

H e {t) = \l \{eV-iv)u s (t,x)\ 2 dx+l [ (p £ (t,x)- p(t,x)) 2 d Xl 

where we have set p e — \u e \ 2 . We find that the time derivative of this modulated 
energy functional is: 



±H e (t) =j J V(divw) -^P e -Y.j d ^ Re (( £a 

j,k 



ivj)u £ (ed k - iv k )u l 



{p e — p) 2 divw. 
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The last two terms are estimated by || X7v(t)\\L^H £ (t). For the first term, write 
e 2 J div (Vt>) ■ V|u £ | 2 = e / div (Vv) • (u £ eVu £ + u £ eVu £ ) 

iv (Vw) • (V(eV - it>)u e + u e (eV - w)u £ ) . 
Since \\u £ {t)\\ L 2 = \\a a \\ L 2 and v € L°°([0, T];W 2, °°), Young's inequality yields: 

l H £ (t)^C(H £ (t)+e 2 ), 

so long as v remains smooth, that is, before shocks appear in (|3.2| . We conclude 
thanks to Gronwall lemma: 

Theorem 3.2. Let n 1, and assume that k = ; V = one? /(y) = V- Assume 
that there exists a smooth function ao independent of e such that 

a% -> a m H s (R n ), Vs > 0. 

Assume also that <f> G C°°(IR n ;M) is suc/i that V0 O G i? s (M") /or a// s ^ 0. T/ien 
i/iere exists T > independent of s > smc/i ttai (|l.ip - (|1.2p feas a unique solution 
u £ G C°°([0,T] x R") n CQO,! 1 ];^) /or a// s > n/2. In addition, ase^O, 

W(zV -iv)u £ \\ 2 Lx{l0 T] . L2) + \\\u e \ 2 - p\\ 2 Loo{[0 T] . L2) =o(e 2 + |||ag| 2 - |ao| 2 ||' 2 ) ■ 

In the above theorem, we have not tried to compute the lowest possible value for 
the Sobolev regularity s given by the proof, nor shall we try in the other sections. 

Remark 3.3. In the more general case where the nonlinearity is f(y) — y a , with 
a G N, a generalization of the above modulated energy functional was introduced in 
3J. In particular, the analogue of Theorem l3.2l is proved. This includes for instance 
the case of the quintic, defocusing nonlinearity. 

One might be afraid that the above result is somehow contradictory with Re- 
mark [3Tj or with the results of [13j . A typical example in |13j . under the assump- 
tions of Theorem l3.21 consists in choosing Oq = ao independent of s, and considering 
v £ solving (TTT ]) -([L2 l) with S o = i 1 + ^ a ) a o {0 < a < 1). Then for t £ of order e a , 

liminf \\u £ {t £ ) -v £ (t £ )\\ L2 > 0. 



Yet, there is no contradiction with Theorem 13.21 the instability mechanism in [13] 
is the appearance of an extra oscillatory factor in v £ . This oscillation shows up 
essentially through a multiplicative factor of the form e l9 ^' x ^ £ . It does not affects 
the modulus of the wave function, and vanishes in the limit e — > of e\?v £ . 

We can therefore conclude that the modulated energy functional shares several 
features with the Wigner measure. It is a rather general tool: in [33], the authors 
consider a nonlinear Schrodinger equation with a boundary condition, aspect which 
apparently cannot be recovered with the approach of E. Grenier recalled in the next 
paragraph. On the other hand, by definition, it ignores the oscillatory phenomena 
that occur at a scale of order s a for < a < 1 (for instance). The next section 
shows how to get a more precise description, under similar assumptions. 



10 



R. CARLES 



3.2. Point-wise asymptotics without potential. In this paragraph, we keep 
assuming V = 0. Note that in (|3.3|) . the supposedly small term on the right hand 
side is of order e 2 , while s should be enough to neglect a term in the limit s — > 0. 
We have seen in i jl.ll that the equation for the phase is obtained after simplification 
by the leading order amplitude. This explains the singular factor on the right hand 
side of (|3.3j) . The main technical ingredient in [32] consists in shifting the source 
term in (|3.3p to the next order, that is, the equation for the amplitude: we now seek 
u e = a 6 e J * l e , where the amplitude a e is complex-valued (even if <Zq is real- valued) , 
3> e is real- valued, and: 



(3.4) 



d t ^ + i|v<r| 2 + /(K| 2 ) = o ; <&f t=0 = 

d t a £ + V$ £ • Va e + -a e A<$> £ = i-Aa e ; af 
2 2 ' !' 



Another originality of this approach lies in the fact that the phase $ e depends on 
e, through the coupling of the two equations. The idea of E. Grenier consists in 
somehow performing the usual WKB analysis "the other way round": first, solve 
(|3.4|) . then show that $ e and a s have asymptotic expansions as e — > 0. In particular, 
this resolves the stability issue pointed out at the beginning of £j3j 

To solve (13.41). consider the new unknown 



u 




n+2 



The system Q3.4p is equivalent to a quasi-linear equation of the form: 



-A ... 

(3.5) d t u e + ^2Aj(u e )djU e = ^Lu £ , with L= I A ... 



e 
2 J 

i=i " \ nxn 



/ v-£ iRea*^ 

and A(u,0=^^i(u)0= v£ |lma^ 

j=l \ 2/' Flea £ 2/'Ima$ w • £T„ 

where /' stands for f'(\a\ 2 ). The system p.5p is hyperbolic symmetric when /' > 0, 
and we can consider the following symmetrizer: 

h 



s 



n i 



which is symmetric and positive for /' > 0. 

Remark 3.4. The argument of /' is morally bounded (this will result from the 
analysis), but may have zeroes: the assumption /' cannot be considered by 
this approach. For instance, justifying a WKB analysis for the quintic, defocusing 
NLS remains an open problem. 

An advantage for this choice of S is that SL remains a skew-symmetric operator: 
the possible loss of derivative caused by the second order operator L does not affect 
the usual energy estimates in H s (J& n ). One can then prove existence and uniqueness 
for ()3.5p in Sobolev spaces of sufficiently large order. Since the last n components 
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define initially, and remain, an irrotational function, this implies that we can solve 
(|3.4|) . The natural limit is given by: 



<9 t $ + - |V$ 



(3.6) 



f(\< 







1 



d t a + V$ • Va + -aA$ = 



*|t=o 



a\t=o — «o 



Local existence in Sobolev spaces for (|3.6I) follows from the same arguments, and 
one has: 

Theorem 3.5. Let n ^ 1, and assume that k — 0, V = and / G C°°(R+; R) wi/i 
/' > 0. Assume that there exists a smooth function ao independent of e such that 



oq in H s 



Vs > 0. 



Assume also that <f) E C°°(R n ;R) is such that Vfo G H s (R n ) for all s^O. Then 
there exists T > independent of e > swc/i that (|l.ip - (|1.2|) /ias a unique solution 
u e = a £ e l ^' E in C°°([0,T] x R") n C([0, T];H S ) for all s > n/2 + 2. Moreover, a £ 
and $ £ are bounded in L°°([0, T]; H s ), uniformly in e G]0, 1] and, /or a/Z s > n/2+X, 
there exists C s such that 

||V$ £ - V$|| ioo([0 T];ff3) + \\a £ - a|| i o O([0iT];ff3) < C s (e + ||ag - a || ff .) . 
Therefore, ||$ e (i) - *(t)|| H . C a t(e+ ||ag - a || H3 ), Vt G [0,T]. 



Theorem 13.51 docs not suffice to describe the asymptotic behavior of u e on the 
time interval [0, T] though: 



ae 1 */ 6 = a e e^ e l e 



i4>/e _ 



(a e -a)e I<r / £ + a( 



Therefore, we have 



ae 



^ |a e - a| + 2|a| 



2e 



Taking the L 2 norm, we infer: 

u £ (t) - a(£)e l * (t)/£ _ < ||a £ (t) - a{t)\\ L 2 + 2||a(t)|| i 2 

(e+||ag-ao|| H .) + M , 



L 2 



$ £ (t) - $(t) 

2e 



< 



for s > n/2 + 1. Even if a 5 , — ao = ©(e^) for AT large, the above estimate shows that 
ae i$/e j g a g 00c j approximation of u e as t — > 0, but not necessarily at time t = T 
for instance. To have a better error estimate, it is necessary to compute the next 
term in the asymptotic expansion of (<p s ,a E ) in powers of e. Assume furthermore 
that there exists a\ G P\ s ^qH s such that 

(3.7) a e = a Q + ea-i + o(e) in Vs > 0. 

For times of order 0(1), the initial corrector ai must be taken into account: 
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Proposition 3.6. Define (a' 1 ', by 

' dM 1] + V$ • V$ (1) + 2/' (|a| 2 ) Rc (aa (1) ) = 0, 

< d t a (1) + V$ • Va (1) + V$ (1) ■ Va + ^a (1) A$ + ^aA$ (1) = ^Aa, 

cbW -n • n (1 1 -n 
L <1> |t=o — > a |t=o^ ai - 

Then aW,$W G L°°([0, T];H S ) for every s ^ ; and 

\\a £ -a- £a (1) |U~( [0 ,T];^) + - $ - e$ (1) ||i-([o,T];^) < C s e 2 , Vs ^ . 

Despite the notations, it seems unadapted to consider $W as being part of the 
phase. Indeed, we infer from Proposition 13.61 that 

u e ~ae l * {1) e l * /e = 0(e). 

L^([0,T];L 2 nL°°) 

Relating this information to the WKB methods presented in £j 1 . 1 1 we would have: 

ao = ae 

Since depends on a\ while a does not, we retrieve the fact that in super-critical 
regimes, the leading order amplitude in WKB methods depends on the initial first 
corrector a\. 

Remark 3.7. The term e l * ' does not appear in the Wigner measure of ae l * ( ' e l */ e . 
Thus, from the point of view of Wigner measures, the asymptotic behavior of the 
exact solution is described by the Euler-type system l|3.2p . 

Remark 3.8. If we assume that ao is real-valued, then so is a. If moreover a\ is 
purely imaginary (for instance, if a\ = 0), then we see that is purely imaginary, 
hence, <J>W = 0. 

Remark 3.9. The proof of Theorem ll.il follows . Consider initial data of the form 
u Q (x) = A-*+ s a (£) , A^O. 

Set e — \ Sc ~ s : e and A go simultaneously to zero, by assumption. Define 
tp e (t,x) = \%- s u(\% +1 - s t,\x) . 

It solves: 

(3.8) ied t r + £ -&r = m 2 r ■, ^ t=0 = a ( X ). 

The idea of the proof is that for times of order 0(1), ip e has become e-oscillatory. 
This is rather clear from (|3.6p : even though <I>| t= o = 0, we have d t &\t=o ^ 0, and 
rapid oscillations at scale e appear instantly. Back to u, this yields the theorem 
(up to replacing ao by | log A| _1 ao). 

To conclude this paragraph, we point out an open problem concerning the time 
T c when shocks appear for ()3.2|) . First, the break-up for (|3.2p does not allow us to 
deduce anything concerning the behavior of the solution of (|3.4p . More generally, 
the notion of caustic in this case is not so clear. Geometrically, as t — > T c , the 
rays for (|3.6p tend to form an envelope. In the linear case / = 0, this geometrical 
phenomenon goes along with an analytical one: 

liminf ||u £ (f)||ioc — > +oo as t — > T c . 

6^0 
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For instance, ||u e (i)||L=° ~ (e+ \T C — 1\) n l 2 for all t in the case of a focal point (all 
the rays meet at one point as t — -> T c ). 

It is not clear at all that a similar phenomenon occurs for (jl.ip when k = 0. 
Suppose for instance that the nonlinearity is cubic, defocusing, f(y) = y, and that 
the initial profile a £ does not depend on e, Oq = a$. The standard conservations of 
mass and energy for nonlinear Schrodingcr equations yield: 

\\u'{t)\\ L , = ||a |U2 = 0(1), 

||eV U £ (i)||2 2 + || U £ (i)||l 4 = ||eVao + ia VMh + IKIll* = ^C 1 )- 

In space dimension n ^ 3, the solution it e remains in ff 1 (R™) for all time, there- 
fore we know that the L 2 and L 4 norms of u £ (t, •) remain bounded by a constant 
independent of e. This suggests that the L°° norm of u e (t, •) may remain bounded, 
if we can somehow inverse the Holder inequality 

iKwiii^Kwiil.iKwiii-. 

One could then distinguish two notions of caustic: a geometrical one (present in all 
the cases), and an analytical one (possibly absent in the highly nonlinear case). 

3.3. Point-wise asymptotics with an external potential. Physical motiva- 
tions may lead to the study of (|l.ip - (|1.2|) when the external potential V is not 
zero. Mathematically, a special role is played by sub-quadratic potentials, as we 
have noticed in §21 see Remark 12.41 We therefore suppose that Assumption 12.11 is 
satisfied. 

The analysis presented in i j3.ll suggests that in this case, we have to consider 
solutions to a compressible Euler equation with (possibly) unbounded external force 
and initial velocity: 

j dtv + v Vv + VV + f'(p)Vp = ; u| t=0 = V0 O , 
(3.9) s „ 

[ dtp + v ■ Vp + pdivv = ; p\ t =o = \ a o\ ■ 

The existence of such solutions is not standard. The naive approach presented in 
[14] consists in resuming the idea of E. Grenier, writing the unknown phase <I> e as 

and considering (|3.4|) where V<jf has replaced V$ £ as an unknown function. This 
procedure is similar to linearizing (|3.4|) in $ £ , around </> e ik- Of course, extra terms 
appear at this stage. Note that the space where we seek $ e is of mixed type: $ £ 
is the sum of a smooth, sub-quadratic (and possibly unbounded) function, and the 
phase 4> e (t, ■) which is sought in Sobolev spaces H s (M. n ). Nevertheless, 4> e must not 
be considered as small, as shown by the analysis of §3.21 

The good news is that the extra terms that have appeared can be treated as 
semi-linear perturbations in the energy estimates. This is due to the fact that 
the phase (/> e ik is sub-quadratic in space. Therefore, the analysis of §3.21 is easily 
adapted: provided that we assume /' > 0, an analogue of Theorem 13. 51 is available. 
Note that unless /' = Const, (in which case the symmetrizer S is constant), we 
need the extra decay assumption on the initial profile: 

xa £ Ds^qH 3 , and xa e -> xa in iP(M"), Vs > 0. 

In particular, a local solution to (|3.9j) is constructed. We refer to [14] for precise 
statements in this case. 
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Remark 3.10. For Schrodinger-Poisson equations in space dimension n 3, the 
idea of E. Grenier was adapted in pjj, under more general geometrical assumptions. 
For instance, solutions that do not necessarily have a zero limit at spatial infinity 
are considered. Under the assumptions of [59] , a point-wise asymptotics of the wave 
function is given, which is more precise that the results in [59] . 



3.4. The case of focusing nonlinearities. Note that in W3.21 the study of 
involves a quasi-linear system whose principal part writes: 



r = a t 2 -div(/'(|^| 2 )v-) 



□ 

This has the same form as the principal part for (|3.6[) , which is the limiting system 
expected in general, whichever formal approach is followed. When /' > 0, we face 
a quasi-linear wave equation. We have pointed out some open problems under the 
weaker assumption /' (a case where loss of hyperbolicity may occur). When 
/' < 0, the above operator becomes elliptic: it does not seem adapted to work in 
Sobolev spaces any more. On the other hand, data and solutions with analytic 
regularity seem appropriate. 

In [30] , P. Gerard works with the analytic regularity, when the space variable x 
belongs to the torus T™, without external potential (V = 0). Note that the only 
assumption needed on the nonlinearity / is analyticity near the range of |ao| 2 . This 
includes the focusing case /' < 0, as well as the defocusing quintic case f(y) = y 2 
for instance. 

The initial phase (f>o is supposed real analytic, and the initial amplitude is analytic 
in the sense of J. Sjostrand [50] : there exist £>Q, A>0, B>0 such that, for all 
j ^ 0, aj is holomorphic in {| lmx\ < £}, and 

\ a j( x )\ < AB j j\ 

Denoting a(t, x) the complex conjugate of a(t, x), P. Gerard considers the system: 
d t v £ = -V s ■ Vv £ - V/ (a ao) , 

1 p j n e 

d t a £ = -v e - Va £ - -a £ divv £ +i-Aa £ (/ {a £ a £ ) - f (a a )) . 

2 2 e 

A solution of the form 

u £ = a £ e^ /£ , a £ (t,x) =J2s j a^(t,x), 

where the sum is defined in the sense of J. Sjostrand, is thus obtained. Setting 

E 

j<l/(C e) 

for Co sufficiently large, the approximate solution v £ satisfies: 

2 

ied t v £ + £ —Av £ = f (\v £ \ 2 ) v £ + O (e~ S/£ ) , 

for some 5 > 0. Essentially, this source term is sufficiently small to overcome the 
difficulty pointed out at the beginning of fj3] for small time independent of e, the 
exponential growth provided by Gronwall lemma is more than compensated by the 
term e~ 8 l £ . We refer to [30] for precise statements and elements of proof. 



v £ = e^ /£ V e 3 a {3) 



SEMI-CLASSICAL LIMIT FOR NLS 



15 



3.5. The integrable case. In the one-dimensional case, n = 1, for a cubic non- 
linearity (f(y) = ±y), the Schrodinger equation is completely integrable. This 
property remains with a time-independent external potential which is a polynomial 
of degree at most two Q] P- 375]. 

In the absence of potential, several papers have studied the semi-classical limit 
for (jl.ip - (|1.2p for the cubic NLS in space dimension one. See for instance [33] in 
the defocusing case, and [35J H3] m the focusing case. A very interesting aspect 
of this approach is that it yields a description of the solution u £ even after shocks 
have appeared for the limiting Euler equation (|3.2p . This description involves theta 
functions, and the so-called Whitham equations (see [55]). In particular, this ap- 
proach seems to confirm the formal discussion of the end of §3.21 in the defocusing 
case, the L°° norm of the solution u £ remains bounded as e — > 0, for all time. 

Unfortunately, it seems that all the results in the integrable case have been 
written in a way that makes any comparison with the other results mentioned 
above very difficult. The last step of inverse scattering is not always performed, 
which should yield a point- wise asymptotics of the wave function u £ . Moreover, 
the spaces in which it would be available are not completely clear. The space 
L^ C (M. X ) seems the most natural candidate. A bridge between the approaches of 
§3.11 and §3.21 on the one hand, and the approaches in the integrable case on the 
other hand, would certainly be welcome in the community of semi-classical analysis 
for nonlinear Schrodinger equations. 



4. Propagation of concentrated initial data 



4.1. Defocusing nonlinearity. We now consider (|1.7p with A > 0. By scaling, we 
may assume A = 1. The general heuristic argument is the following. For t close to 
zero, the solution u £ remains concentrated near the point xq, at a scale of order e. 
Since the potential V does not depend on e, we have Vu £ ~ V(xo)u £ : the potential 
can be considered as constant at leading order. Introduce the function ip £ given by 
the scaling 



v. £ {t,x) 



1 



i/2 



r - 



t X-X \ i{x .^ /e-V{x )t/ S ) 



The Cauchy problem (jl.7p is equivalent to: 
id t ip £ + -Aip e = {V(x +ex)~V{x ))ij £ + |V £ | 2< > e 5 ifj £ {0,x) = R(x)e tx<0 . 

xp, where ip is independent of e 



The above argument suggests that we have ip' 
and solves: 



(4.1) 



ip(0,x) = R(x)e tx - io . 



Under suitable assumptions on a and R, there is scattering for this equation (see 
e.g. [13 DI1 S3]): there exist tp± € L 2 (R n ) such that 



(4.2) 



as t — > ±oo. 



The standard asymptotics of the free Schrodinger group e J 2 A then yields: 

e t\x\ 2 /(2t) 



fl>(t, 



t->±oo (it)"/ 2 



-&(?). 
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where the Fourier transform is given by 

Back to u E , this yields, for t ^> e and so long as we consider the external potential 
as constant: 

(4.3) u'(t,x) ~ j^ e iM/ s -v ixoWs) _ 
Indeed, we have the following rigorous result: 

Proposition 4.1 ([16], Proposition 6.3). Let V satisfying Assumption UOJ Let 
JJeS:= H 1 n TiH 1 ), and 

2-n + \Ai 2 + 12n + 4 2 

^ a < 

An n-2 

Then for any A > 0, the following holds: 

1. There exists e(A) > such that for < e < e(A), the initial value problem (|1.7|) 
has a unique solution u 6 € C([— Ae, Ae]; S). 

2. This solution satisfies the following asymptotics, 

limsup sup \\u £ (t) — v £ (i)|| L 2 = , 

£^0 |t|<Ae 

where v e is given by v e (t,x) — n ^ 2 ip 

and tp S C(R; S) is given by (|4.ip . 

A transition is expected to occur in the above boundary layer, that is for \t\ = Ae 
and A ^> 1. The heuristic argument consists in saying that because of dispersion for 
ip, the external potential V can no longer be considered as constant. On the other 
hand, and for the same reason, the nonlinearity ceases to be relevant at leading 
order: for e <C ±i ^ T, we expect u 6 ~ u§_, where 

(4.4) ied t u% + — Aui = Vu% ; u e (0, x) = — ^?/>± (-) e ^o<o/e-v(x )t/e) ^ 

and tp± are given by (|4.2[) . The value of T is not arbitrary: the asymptotic behavior 
of u|_ involves the classical trajectories associated to V. These trajectories may 
refocus at one point; this is the case when V is an isotropic harmonic potential for 
instance. 

Proving the above asymptotics for e <C ±t ^ T is actually an open problem for 
general potentials satisfying Assumption 12.11 even for time- independent potentials. 
It has been proved when V = V(x) is exactly a polynomial of degree at most two, 
in [TT] for the case of refocusing(s), and in [16] for the complementary case. 

The restriction to this class of polynomial potentials is certainly purely technical, 
and we know explain it. The proof of the asymptotics for e <C ±t ^ T relies on the 
use of operators well suited to the propagation of classical trajectories associated 
to V. In the linear setting, good candidates to meet this requirement are given by 
the action of Heisenberg derivatives (see e.g. [48]): 

U e {t)eVU £ (-t) and U s (i)?—^-U e (-t) , where U e {t) = e _i *(r^ A+y ). 
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The main technical remark in I16j is that when V is a polynomial of degree 
at most two, then the above two Heisenberg derivatives are very interesting for 
nonlinear problems too. Indeed, we can find p = p(t), and <fi = ef>(t,x) real- valued, 
such that, for instance: 

(4.5) u £ (t)^—-^-u e {-t) = P {ty t ' < - t ' x)/£ v (V^*)/^ . 

In [16] , it is proved that an operator of the form of the right hand side of (|4.5p 
commutes with U E (t) if and only if V is a polynomial of degree at most two, and <f> 
solves the eikonal equation (|2.ip . The fact that an Heisenberg derivative commutes 
with the group U £ (t) is a straightforward consequence of its definition. The right 
hand side of (14. 5|) implies two important things: 

• This Heisenberg derivative acts on gauge invariant nonlinearities G(|u| 2 )u 
like a derivative. 

• Weighted Gagliardo-Nirenberg inequalities are available, of the form 



L 2 

To illustrate the use of these properties, we recall [TTJ Corollary 1.3]: 
Proposition 4.2. Let fieE. Assume that u e solves (|1.7p with x = £ = and 

Let i/)± = W± l R be given by (|4.2[) (upon suitable assumptions on a, see e.g. 
Prop, \4-l\l - Then for any 2 < r < -^tt, the following asymptotics holds in L 2 n L r : 

/ ^ \ n/2 2 

• IfO < t < 7T. then u e (t,x) ~ [ ) ib+ (— e l ^n . 

v ' e-o \isint J \smtJ 

• If -tt <t < 0, thenvL e (t,x) ~ f— !— ^) ^sifil*. 

1 y e^o \isint J v Vsini/ 

Remark 4.3. The result of [5] shows that in the above case, Wigner measure is not a 
good tool to characterize the behavior of u e . More precisely, we can find i?i , i?2 £ E 
such that the Wigner measures for the corresponding solutions uf and u| coincide 
at time t = — tt/2, but are different at time t = ir/2. The crossing of a focal point 
may lead to an ill-posed Cauchy problem as far as Wigner measures are concerned. 



We see that the formal asymptotics (|4.3|) is valid only in the transition regime 
t = As, with A > 1. For larger times, the trigonometric functions in the above 
result account for the dynamical influence of the harmonic potential. 

In the above case of an isotropic harmonic potential, the above result can be 
iterated in time. Recall that the (nonlinear) scattering operator S associated to 
(|4.1[) maps tp_ to tp + , given by (|4.2|) . 



Corollary 4.4. Under the assumptions of Proposition pTz} consider k € N. For 
kir < t < (k + and 2 < r < -^z^, the following asymptotics holds in L? n L r : 



u £ (i,z) — — F -S k i!+ — e 1 ^ 

v e-o |sint|"/2 v+ Vsint/ 



where S denotes the k th iterate of the scattering operator S . 
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The phase shift e~ mk % corresponds to successive Maslov indices: this is a linear 
phenomenon [25) . On the other hand, we see that a nonlinear phenomenon occurs 
at leading order at time t = kit, which is measured by the scattering operator S. 

4.2. Focusing nonlinearity. When A < in (|l.ip - (|1.3p . we assume similarly that 
A = — 1. We let R = Q, the unique positive, radially symmetric solution of: 

~AQ + Q= \Q\ 2 °Q. 

Now, the focusing nonlinearity is an obstruction to dispersive phenomena. The 
solution u e is expected to keep the ground state as a leading order profile. Never- 
theless, the point where it is centered in the phase space, initially (i£o,£o); should 
evolve according to the Hamiltonian flow (|2.2[) . In the absence of external potential, 
V = 0, we have explicitly: 



U s (t,x) = Q ( X e ix<(t)/e+iB(t)/e 



where (x(t),£(t)) = (x +f£oi£o) solves (|2.2I) with initial data (xq,£q), and 9(t) = 
t— £|£o| 2 /2. When V is not trivial, seek u e of the form of a rescaled WKB expansion: 

i<f>(t,x)/e 




u s (t,x)~ }eW 3 -, ^ e 



Note that this scaling meets the exact result of the case V = 0. Plugging this 
expansion into (|1.1|) - (|1.3|) and canceling the O(e ) term, we get: 

id t Uo + ^AU + U (-d t <t> - i|V0| 2 -V+ \U \ 2a ^j - i (x(t) - V</>) ■ VU = 0. 

Impose the leading order profile to be the standing wave given by 

U {t,x) = e it Q(x). 
Then the above equation becomes: 

U (-d t cf> - i|V0| 2 -v\-i (x{t) - • V(7 = 0. 

Since Uoe~ lt is real- valued, and since we seek a real-valued phase (f>, this yields: 
d t( f>+^\V^\ 2 + V = ; cl>(0,x)=x-Z . 
x{t) = V<t>(t,x). 

The first equation is the eikonal equation (|2.1[) . We infer that we have exactly 

(*,*(<))=£(*). 

The form of Uo and the exponential decay of Q show that we can formally assume 
that x = x(t) + 0(e). In this case, 

V<0(t, a:) = V<0 (t, x(t)) + O(e) = f (t) + O(e) = x(t) + O(e). 

Thus, we have canceled the 0(e°) term, up to adding extra terms of order e, that 
would be considered in the next step of the analysis, which we stop here. Back to 
u E , this formal computation yields 

u ^t,x) ~q( ^1^1) e W,») „q( „to-6(t)/ B -Hfl(t)/ e 
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where 6{t) = t (l - |£o| 2 / 2 - V(x )) + / x(s) ■ VV(x(s))ds. 

Jo 

To give the above formal analysis a rigorous justification, the following assump- 
tions are made in [35]: 

Assumption 4.5. The nonlinearity is L 2 -subcritical: a < 2/n. 

The potential V = V{x) is real-valued, and can be written as V = Vi + Vi, where 

• Vie w 3 >°°(w i ). 



d a Vn £ W 2,00 (R™) for every multi-index a with \a\ = 2. 



For instance, V can be an harmonic potential. 

Theorem 4.6 ([38]). Let xq,^o € K™. Under Assumption \4-5\ the solution u e to 
(|l.l|) - (|1.3|l with R~Q can be approximated as follows: 



* M - Q (^®.) ,-«.»/.-«.>, +0(£) hIR 



Ixe = ^ll/lli3 + ^l|v/||i 2 . 



where (x(t),^(t)) is given by the Hamiltonian flow, the real-valued function 9 e de- 
pends on t only, and X £ is defined by the norm 

Remark 4.7. The assumption a < 2/n is crucial for the above result to hold. Indeed, 
if a = 2/n and V is the isotropic harmonic potential 

then we have explicitly, when xq = £o = (see [101 [38]): 

u e (t,x) = - l — T rQ(—^-)e^-^^ t , 0<t<£, 

v ' ; (cost) n /2^ Vecost/ 2 

so the profile Q is modulated as time evolves, in a fashion similar to §4.11 

The proof of the above result heavily relies on the orbital stability of the ground 
state, which holds when a < 2/n. For v £ iJ 1 (R"), denote 

m = l\\Vv\\»-^\M%t& 

The ground state Q is the unique solution, up to translation and rotation, to the 
minimization problem: 

£(Q) = w£{£{v) ; v E if 1 (M n ) and ||«|| £a = \\Q\\^}. 

The orbital stability is given by the following result: 

Proposition 4.8 ([54]). Let a < 2/n. There exist C, h > such that if 4> e H 1 (R n ) 
is such that \\^>\\l 3 = ||Q||l 2 an d ~~ £{Q) < h, then: 

inf U-e ie Q(--y)\\ 2 H1 < C(£(0) - . 

The strategy in [38] consists in applying the above result to the function 

v e (t, x) = u £ (t, ex + x(t)) e -^*+x(t)K(i)/ e _ 
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For A > sufficiently large, let \ be a smooth non-negative cut-off function, sup- 
ported in {x £ M n ; |a;| ^ 2A}, and constant equal to 1 in {x € K n ; |x| ^ A}. 
Introduce the error estimate rj £ (t) given by r] s = r\\ + t}\ + V3 + Vl j where: 

»7i(*)= / a;x(x)TO £ (t,a;)da;- ||Q||| 2 a;(i), 

JR™ 
JR™ 

7?i(t)= / f(*,^-iiQii| 2 ew, 

JR™ 
JR" 

m £ {t,x) = -^\u e (t,x)\ 2 ; £ e (t,aO - -L Im (H e Vii £ ) . 

Noting that ?f (0) = 0{e 2 ), the proof in [38] shows that rf (t) = C(e 2 ) for f e [0,T ] 
for some Tq > independent of e. The proof eventually relies on Gronwall lemma 
and a continuity argument. In order to invoke these arguments, S. Keraani uses 
Proposition ^. 81 and the scheme of the proof of J. Bronski and R. Jerrard [6], based 
on duality arguments and estimates on measures. Finally, the time To given by the 
proof depends only on constants of the motion, so the argument can be repeated 
indefinitely, to get the Lf£ c estimate of Theorem 14.61 

In the particular case where the external potential V is an harmonic potential 
(isotropic or anisotropic), the proof can be simplified. We invite the reader to pay 
attention to the short note [37], where this simplification is available. 

The phase shift 8 £ in Theorem 14.61 is not known in general. It is easy to guess 
from the arguments given above that in the proof given by S. Keraani, it stems from 
the use of Proposition 14.81 On the other hand, as noted in [38 , a time-dependent 
phase shift does not alter the Wigner measure of u 6 , which is an important physical 
quantity. 

Acknowledgments. The author is grateful to Thomas Alazard for his careful 
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